207 lines (207 with data), 5.8 kB

**
**BONE GEOMETRY FILES
*include, input=bone_femur.inp
*include, input=bone_tibia_fibula.inp
**
****BONE PROPERTY FILES, including masses, rigid body reference nodes****
*include, input=props_bone.inp
**
****CARTILAGE GEOMETRY FILES****
*include, input=cart_femur.inp
*include, input=cart_tibia.inp
**
****CARTILAGE PROPERTY FILES****
*include, input=props_cart.inp
**
****GS TF LOAD CELL FILE****
*include, input=loadcell_tf.inp
**
**************************************************************
**
****LIGAMENT PARAMETERS****
****COMMENT OUT SPECIFIC LIGAMENTS TO MODEL DIFFERENT RESECTION STAGES FROM IN VITRO TESTS****
*include, input=lig_parameters.inp
**
*include, input=lig_als.inp
*include, input=lig_acl.inp
*include, input=lig_lcl_united.inp
*include, input=lig_mcl.inp
*include, input=lig_pcap.inp
*include, input=lig_pcl.inp
*include, input=lig_pfl_slipring.inp
**
**
****BEHAVIORAL PROPERTIES OF LIGAMENTS****
*include, input=props_als.inp
*include, input=props_acl.inp
*include, input=props_lcl_united.inp
*include, input=props_mcl.inp
*include, input=props_pcap.inp
*include, input=props_pcl.inp
*include, input=props_pfl_slipring.inp
**
****LIGAMENT ATTACHMENTS****
*include, input=attach_ligs_slipring.inp
**
*include, input=dummy.inp
**
**************************************************************
****GROOD SUNTAY NODES AND AXIS DEFINITIONS****
*include, input=gs_nodes.inp
*include, input=gs_kin_axes.inp
**
**************************************************************
**
**** INTERACTION DEFINITIONS ****
*SURFACE INTERACTION, NAME=RIGID_INTERACT
*SURFACE BEHAVIOR,PRESSURE-OVERCLOSURE=LINEAR
3.2
**
**
**************************************************************
**** LIGAMENT ATTACHMENT SITE TRANSLATION ****
*include, input=lig_attach_translation.inp
**
**************************************************************
**** STEP 1: SETTLE MODEL INTO POSITION AT DESIRED FLEXION ANGLE ****
*step
*dynamic, explicit, direct user control
2.0e-04, 1.0
**
**** TIME INTERVAL CAN BE INCREASED FOR THE SETTLING STEP TO REDUCE SIZE OF OUTPUT FILES
*OUTPUT,FIELD,VARIABLE=PRESELECT, TIME INTERVAL=0.2
**
*INCLUDE, INPUT=contact_rigid_CP.inp
**
*include, input=amps_0settling_step.inp
**
**** SOME BOUNDARY CONDITIONS: FIX THE FEMUR RIGID BODY REACTION NODE, 
**** FIX THE MATERIAL FLOW DEGREE OF FREEDOM OF NODES INCLUDED IN THE PFL SLIPRING
*boundary, op=new
1100000, encastre
8009, 10,10, 0.0
8010, 10,10, 0.0
8011, 10,10, 0.0
8001, 10,10, 0.0
8002, 10,10, 0.0
8005, 10,10, 0.0
**
**
**** MOVE JOINT TO DESIRED FLEXION ANGLE ****
*connector motion, amplitude=flexion
axis_TF_ML, 4, 0.0174
**
**** INITIAL VAL-VAR POSITION (DEG) ****
*connector motion, amplitude=VV_initial
axis_TF_AP, 4, 0.0174
**
**** INITIAL IE POSITION (DEG) ****
*connector motion, amplitude=IE_initial
axis_TF_SI, 4, 0.0174
**
**** INITIAL ML POSITION ****
*connector motion, amplitude=ML_initial
axis_TF_ML, 1, 1
**
**** INITIAL AP POSITION ****
*connector motion, amplitude=AP_initial
axis_TF_AP, 1, -1
**
**** INITIAL SI POSITION - CURRENTLY NOT USED ****
***connector motion, amplitude=SI_initial
**axis_TF_SI, 1, 1
**
** ADD A SMALL SI LOAD ON FEMUR ****
*connector load, amplitude=load_SI
axis_TF_SI, 1, -1.0
**
**
**** FIX THE LOADCELL IN ITS POSITION WITHIN THE TIBIA SO THAT IT CAN MEASURE LOADS ****
*connector motion
loadcell_TF, 1, 0.0
loadcell_TF, 2, 0.0
loadcell_TF, 3, 0.0
loadcell_TF, 4, 0.0
loadcell_TF, 5, 0.0
loadcell_TF, 6, 0.0
**
**** OUTPUT FILES FROM STEP 1
*include, input=output_step1_kinematics.inp
*include, input=output_step1_ligs_cruciates_slipring.inp
*include, input=output_step1_loadcell_tf.inp
*include, input=output_step1_cartilage.inp
**
**
*end step
**
**************************************************************
**** STEP 2: AT DESIRED FLEXION ANGLE, MATCH THE EXCURSION PROFILES ****
*step
*dynamic, explicit, direct user control
2.0e-04, 0.6
**
*OUTPUT,FIELD,VARIABLE=PRESELECT, TIME INTERVAL=0.1
**
*include, input=amps_0flex_anterior.inp
**
**** SOME BOUNDARY CONDITIONS: FIX THE FEMUR RIGID BODY REACTION NODE, 
**** FIX THE MATERIAL FLOW DEGREE OF FREEDOM OF NODES INLCUDED IN THE PFL SLIPRING
*boundary, op=new
1100000, encastre
8009, 10,10, 0.0
8010, 10,10, 0.0
8011, 10,10, 0.0
8001, 10,10, 0.0
8002, 10,10, 0.0
8005, 10,10, 0.0
**
**** Designate the axes about or along which rotation or translation should occur ****
*connector load, op=new, amplitude=maintain_load_SI
axis_TF_SI, 1, -1.0
**
*connector motion, op=new, amplitude=maintain_flexion
axis_TF_ML, 4, 0.0174
**
***connector load, op=new, amplitude=TF_VVLoad
**axis_TF_AP, 4, 1
**
***connector load, op=new, amplitude=TF_IELoad
**axis_TF_SI, 4, 1
**
*connector load, op=new, amplitude=TF_APLoad
axis_TF_AP, 1, -1
**
***connector motion, op=new, amplitude=TF_VV
**axis_TF_AP, 4, 0.0174
**
***connector motion, op=new, amplitude=TF_IE
**axis_TF_SI, 4, 0.0174
**
*connector motion, op=new, amplitude=TF_ML
axis_TF_ML, 1, 1 
**
***connector motion, op=new, amplitude=TF_AP
**axis_TF_AP, 1, -1
**
**** SI KINEMATIC CONTRAINTS COMMENTED OUT TO ALLOW MODEL FREE MOTION ****
***connector motion, op=new, amplitude=TF_SI
**axis_TF_SI, 1, 1
**
**
**** FIX THE LOADCELL IN ITS POSITION WITHIN THE TIBIA SO THAT IT CAN MEASURE LOADS ****
*connector motion, op=new
loadcell_TF, 1, 0.0
loadcell_TF, 2, 0.0
loadcell_TF, 3, 0.0
loadcell_TF, 4, 0.0
loadcell_TF, 5, 0.0
loadcell_TF, 6, 0.0
**
**** OUTPUT FILES FROM STEP 2 ****
**
*include, input=output_step2_kinematics.inp
*include, input=output_step2_ligs_cruciates_slipring.inp
*include, input=output_step2_loadcell_tf.inp
*include, input=output_step2_cartilage.inp
**
*end step